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Despite the fact that the solar neutrino flux is now well-understood in the context of matter- 
affected neutrino mixing, we find that it is not yet possible to set a strong and model-independent 
bound on solar neutrino decays. If neutrinos decay into truly invisible particles, the Earth-Sun 
baseline defines a lifetime limit of r/m > 10 -4 s/eV. However, there are many possibilities which 
must be excluded before such a bound can be established. There is an obvious degeneracy between 
the neutrino lifetime and the mixing parameters. More generally, one must also allow the possibility 
of active daughter neutrinos and/or antineutrinos, which may partially conceal the characteristic 
features of decay. Many of the most exotic possibilities that presently complicate the extraction of a 
decay bound will be removed if the KamLAND reactor antineutrino experiment confirms the large- 
mixing angle solution to the solar neutrino problem and measures the mixing parameters precisely. 
Better experimental and theoretical constraints on the 8 B neutrino flux will also play a key role, 
as will tighter bounds on absolute neutrino masses. Though the lifetime limit set by the solar flux 
is weak, it is still the strongest direct limit on non-radiative neutrino decay. Even so, there is no 
guarantee (by about eight orders of magnitude) that neutrinos from astrophysical sources such as a 
Galactic supernova or distant Active Galactic Nuclei will not decay. 

PACS numbers: 13.35.Hb, 14.60.Pq, 26.65.+t FERMILAB-Pub-02/061-A 



I. INTRODUCTION 

Since solar neutrinos have been detected with roughly 
the expected flux, it appears that they do not decay over 
the 500 s x c distance to Earth. Furthermore, neutrinos 
from SN1987a were also detected in reasonable numbers 
on a much longer baseline of (5 x 10 12 s) x c. Decay 
will deplete the flux of neutrinos of energy E and mass 
m over a distance L by the factor 

eXP ("^)^ eXP H X v)' 

where r is the rest-frame lifetime and we use c = 1 units 
from now on. In Table I, we list representative r/m scales 
for various neutrino sources. 

In this paper, we critically assess what the best limits 
on neutrino decay are. We find that it is not yet possi- 
ble to set model-independent bounds, even for the well- 
measured solar neutrinos. We discuss how decay limits 
can be improved in the future. 

Though in the past neutrino decay was frequently dis- 
cussed in terms of flavor eigenstates, the lifetimes of neu- 
trinos are only well-defined for mass eigenstates (a flavor 
state does not have a definite mass, lifetime, or magnetic 
moment). Since we now know that mixing angles are 
large, this distinction is essential. 

Therefore, in considering the decay of neutrinos from 
the Sun and SN1987a, one has to properly assess the 
mass eigenstate content of the fluxes. The SN1987a data 
can be reasonably explained by saying that the expected 
flux of T>\ made it to Earth and was detected as v e via 
the charged-current reaction v e p — > e + n. The compo- 
nent would only have been detectable in neutral-current 



TABLE I: Representative scales for neutrino lifetimes, taken 
simply as r/m ~ L/E. The top three entries correspond to 
present data, and the lower two to possible future data (after 
neutrinos of all flavors have been observed from a Galactic 
supernova and neutrinos from astrophysical sources like active 
galactic nuclei (AGN) have been observed in km 3 detectors). 



Neutrino source 


L/E r/m 


(s/eV) 


Accelerator 


30 m/10 MeV 


1Q -14 


Atmosphere 


10 4 km/300 MeV 


10 -io 


Sun 


500 s/5 MeV 


10" 4 


Supernova 


10 kpc/10 MeV 


10 5 


AGN 


100 Mpc/1 TeV 


10 4 



reactions, and the SN1987a data are consistent with no 
neutral-current events. If v\ is (as suggestively labeled) 
the lightest mass eigenstate, then it would not be kine- 
matically allowed to decay, so that the decay limit from 
SN1987a would be meaningless. 

Presently, the best explanation of the solar neutrino 
problem is large mixing angle (LMA) Mikheyev-Smirnov- 
Wolfenstein (MSW) transformation of v e to v^,v T . Be- 
sides being the best oscillation-parameter fit, LMA also 
provides a "good" fit to all of the solar neutrino data in 
terms of an acceptable chi-squared. Without the effects 
of oscillations, the solar neutrino flux is only understood 
at the factor of two level. Including the effects of oscil- 
lations, the total flux of all flavors is better understood, 
and the hope is that much more stringent decay limits 
can be derived. A similar approach was used to set the 
strongest direct neutrino magnetic moment limit [Q. 

If LMA is the correct scenario, then solar neutrinos in 



the 5 — 15 MeV energy range are created as nearly pure 
v 2 mass eigenstates. Furthermore, their propagation in 
the Sun is completely adiabatic, so that they emerge as 
pure v 2 eigenstates, where 777,2 > m\. Since the charac- 
teristic signature of decay is its energy dependence, we 
restrict our attention to this energy range, over which 
Super-Kamiokande (SK) Q and the Sudbury Neutrino 
Observatory (SNO) || have measured energy spectra. 
This data may thus be used to search for the signatures 
of the kinematically allowed v 2 — > v x + X decay (the 
decay modes are discussed in detail below). 

While neutrino decay was an early proposed explana- 
tion of the solar neutrino problem ||, |5| , in this paper we 
will take the point of view that the LMA solution specifies 
the correct basic picture and will consider v 2 — * V\ + X 
decay as a perturbation, with the goal of placing a limit 
onr/m (since m is unknown, only the quantity r/m can 
be constrained). Ongoing experiments will soon confirm 
or refute the LMA solution. 

While there are strong limits on radiative neutrino 
decays 0, ffl], in this paper we consider only the so- 
called "invisible" decays, i.e., decays into possibly de- 
tectable neutrinos (or antineutrinos) plus truly invisible 
particles, e.g., light scalar or pseudoscalar bosons. For 
these modes, the lifetime limits are very weak indeed; 
see Table I. There are also limits from cosmology, the 
strongest of which use the cosmic microwave background 
data ||. The effect searched for arises because of the 
transfer of nonrelativistic energy density in the parent 
neutrinos (m > 10 eV) to relativistic energy density in 
the daughter neutrinos (m ~ 0). Neither the large mass 
of the parent neutrino nor the huge mass splitting is sup- 
ported by present data, so these cosmological limits on 
non-radiative neutrino decays are inapplicable. 

While detection of daughter neutrinos or antineutrinos 
has been considered in the literature, it has nearly always 
been with the assumption that mi <C m 2 , and that the 
daughter neutrino carries half the energy of its parent 
in a two-body decay. These assumptions are not gener- 
ally valid. Since oscillation experiments do not determine 
the overall mass scale, the lightest mass eigenvalue mi is 
unknown. However, for fixed mi, the masses 7712 and 
777,3 are determined by the measured mass-squared differ- 
ences. Thus 



m 2 = \J mj - m\ + ml = y m 2 + Sm 2 sol , (1.2) 

where solar-neutrino data give 5m 2 ol ~ 4 x 10~ 5 eV 2 , 
and similarly, 

m 3 = \J ml + Sm 2 sol + Sm 2 atm , (1.3) 

where atmospheric-neutrino data give Sm 2 ^ ~ 3 x 1CP 3 
eV 2 . All three masses are shown in Fig. ffl, illustrating 
that the masses are nearly degenerate unless the overall 
mass scale is tiny. The accepted wisdom is that neutrino 
masses should be strongly hierarchical, e.g., in see-saw 
models, proportional to the squares of charged-lepton 
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FIG. 1: The values of the neutrino mass eigenvalues as a 
function of the unknown smallest mass mi . When mi is spec- 
ified, the measured solar and atmospheric Sm 2 values fix mi 
and 772,3. For an inverted hierarchy (<5mj; tm < 0), the two up- 
per eigenvalues are very nearly degenerate at the position of 
the curve labeled 7713. The present bound on neutrino mass 
from tritium experiments is about 2 eV ||; by the construc- 
tion above, it applies to all three mass eigenvalues. 



masses. Note that the ratio 7713/7712 < 10 for all m\. 
Naively, this argues against a simple see-saw mass pat- 
tern and argues for the idea that the mass ratios can be 
small, and quite possibly degenerate. 

If the neutrino masses are degenerate, then the daugh- 
ter neutrino carries nearly the full energy of the parent 
neutrino in any reference frame (in the rest frame, a v 2 
at rest decays to a v\ nearly at rest). This completely 
alters scenarios in which active daughters are detected. 



II. NEUTRINO DECAY MODELS 

Non-radiative neutrino decay may arise through the 
coupling of the neutrino to a very light or massless parti- 
cle, such as a Majoron ]l0| . Majoron models typically 
have tree-level scalar or pseudoscalar couplings of the 
form 

£ = gijViVjX + hijV i 'j 5 VjX + h.c, (2.1) 

where x is a massless Majoron, which does not carry 
definite lepton number. For the couplings specified by 
Eq. (pTTl) , the decay rates into neutrino and antineutrino 
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daughters are given by Ref. |Uj 
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where x = 7712/7711, and we have dropped the subscripts 
on the coupling constants. The decay widths in this sec- 
tion are defined in the laboratory frame, so the relation 
to the rest-frame lifetimes quoted elsewhere is 
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(2.4) 



In the limit of hierarchical neutrino masses, 7772 3> 7771, 
the case that has received the most attention to date, the 
decay rates are equal: 



(ff 



32tt 3 



^- (2-5) 



In the opposite limit, where 7772 
Sm 2 0), we find instead: 



777 1 (but keeping 
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(2.6) 



The important point to note is that in these Majoron 
models, both neutrino and antineutrino decay products 
may be produced, but the relative weight of the two decay 
modes depends strongly upon the mass hierarchy. 

In the simplest versions of these models, the neutrino 
masses are proportional to these coupling constants and 
hence the neutrinos are exactly stable, as the matrix of 
coupling constants is diagonal in the mass basis. Even in 
this case, the neutrinos may have finite lifetimes in mat- 
ter, as the rotation of the mass basis in matter will lead 
to non-diagonal couplings between matter mass eigen- 
states Jl^, In the most general case, the basic models 
can easily be modified to permit non-diagonal couplings 
in vacuum jl4[ |. 

While there are a huge variety of Majoron models 
in the literature including, for instance, "charged" Ma- 
joron and "vector Majoron" jl6| models, we will not 
restrict our attention to any particular model. For exam- 
ple, instead of a Majoron model per se, we can consider 
the couplings of neutrinos to a very light gauge boson. 
Similarly, we make no assumption about the relationship 
between the neutrino masses and the couplings that give 
rise to decay, or whether the neutrinos are Dirac or Ma- 
jorana particles. 

We thus take a purely phenomenological point of view 
and consider any possible tree-level coupling between the 



neutrinos and a very light or massless particle. We shall 
consider the cases where the decay products are active 
or sterile neutrinos or antineutrinos, plus an "invisible" 
particle. While one could imagine models where such a 
coupling is absent at tree level (as with radiative decay 
in the Standard Model, for example) , models of this sort 
are of less interest, as any coupling arising only at loop 
level is likely to lead to a very small decay rate. 

B ound s on neutrino-Majoron couplings of the form in 
Eq. (2.1) may be obtained from considering their effects 
on the two-body leptonic decays of 7r and K mesons at 
rest. A nonzero coupling allows the final neutrino to also 
appear as a neutrino or antineutrino, plus a Majoron. 
This increases the decay rate and also smears the mo- 
mentum distribution of the charged lepton. The bounds 
obtained are approximately g 2 < 1CP 4 Jl7], ^8), and are 
reasonably model- independent. 1 (Note that in this sec- 
tion we denote by g either a scalar or pseudoscalar cou- 
pling) . A considerably more stringent bound may be de- 
rived from limits on ncutrinoless double beta decay with 
Majoron emission [jl9|. However, the limit of g 2 < 10~ 8 
applies only to the g ee element of the coupling matrix and 
does not directly translate into a bound on the parameter 
of interest, namely gi\. 

Translated into a bound on neutrino lifetimes, the 
meson-decay bound on the coupling of g 2 < 10~ 4 be- 
comes 



— > 3 x 1CT 5 s/eV 



1(T 5 eV 2 
Sm 2 



(2.7) 



where we have used Eq. fl2.5p . With the solar Sm 2 ~ 4 x 
10 -5 eV 2 , the lifetime limit obtained allows substantial 
decay of solar neutrinos (from Table I, the natural scale of 
the problem is r/m ~ 10~ 4 s/cV). If the mass hierarchy 
is inverted, then solar neutrino decays can also occur with 
Sm 2 ~ 10~ 3 eV 2 , in which case the derived lifetime limit 
is 100 times weaker, and hence not useful. 

We emphasize again that it is not our intention to re- 
strict our attention only to the case of Majoron models. 
We discuss these models only as a concrete example of 
the types of decay modes that may be expected. In Sec- 
tion 111 we take a quite general perspective and discuss 
specific decay modes as illustrative, model-independent 
examples. 



III. SOLAR v 2 DECAY 

We consider the decay of neutrinos during their jour- 
ney from the Sun to Earth, and neglect decay inside the 



1 These bounds do not apply to the g TT element of the coupling 
matrix which, given the large neutrino mixing angles, will con- 
tribute to all elements of g in the mass basis. Whether it is likely 
that g TT would be significantly greater than all other (flavor ba- 
sis) elements g is a model-dependent question of naturalness. 
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Sun, since the distances are 500 s and 2 s, respectively. 
Decay rates can be increased in matter Q by the greater 
phase space provided by the matter enhancement to Sm 2 . 
However, for LMA-type solutions, this effect is not large. 
If decay in the Sun were significant it would imply that 
our understanding of the solar neutrino flux is funda- 
mentally flawed; KamLAND and future solar-neutrino 
experiments will check this possibility. 

As noted, since an identifying characteristic of decay 
is its particular energy dependence, we consider only the 
SK and SNO spectral data. These experiments are sen- 
sitive to 5 — f5 McV neutrinos from 8 B beta decay (while 
the normalization uncertainty is ~ 20%, it is common 
to SK and SNO). While one might think that the low- 
est energy neutrinos are the best suited for testing decay 
effects, there are two important caveats. First, the gal- 
lium and chlorine radiochemical experiments only mea- 
sure energy-integrated rates, and receive contributions 
from several different neutrino sources with uncertain 
relative normalizations. Exotic effects on the neutrino 
survival probability can be hidden in this data. Second, 
for the LMA solutions, the v\ to v-i ratio rises at low 
energies, and as noted, only the heavier V2 can decay. 

In the energy range covered by the SK and SNO spec- 
tra, the solar neutrino flux is dominated by the 8 B neutri- 
nos. To a very good approximation, these neutrinos are 
produced at the solar center where the matter potential 



is 



C = V2G F N t 



2E 
bra 2 



15.3 



E \ /10" 5 eV 2 



10 McV 



Sm 2 



(3-1) 



The initial v\ and fractions of the flux are specified by 



P l Vl = cos 2 9 m , 



P' 



where the matter mixing angle given by 

• 2 ofl sin 2 20„ 
sm 26 m — 



sin 2 26 v + (cos26> m - C) 2 



(3.2) 



(3.3) 



As representative values 
we choose 8m 2 



4 x 



and 6 V is the vacuum angle 
in the LMA allowed region |20 

10" 5 eV 2 and sin 2 26 v = 0.9. The MSW transformation 
is always adiabatic in the solar neutrino energy range 
for the LMA solutions, so the fluxes of V\ and v-i are 
unchanged at the solar surface (neglecting the tiny decay 
probability over the solar radius). The final v\ and v-i 
fluxes at Earth, taking decay into account are 



PJ 

pf 

V2 



= P l 

„ f L m 

= P , exp x — 

1/2 y V E T 



(3.4) 



Note that the small variation in the Earth-Sun distance 
cannot be exploited, as its effects are washed out by the 
energy resolution. It is convenient to think of a typical 



energy for decay (that makes the argument of the expo- 
nential unity); in convenient units, this is 



E de 



1 MeV 



5 x 10~ 4 s/eV 
r/m 



(3.5) 



The electron neutrino survival probability at Earth is 
thus 



P{v e -» v B ) 



Pi, cos 2 



0, + Pi 



The muon neutrino appearance probability is 
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P£ sin 2 



P/ 2 cos 2 
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(3.7) 



Because of decay, P(v e — > v e ) + P(v e — ► v^) < 1 (at 
this point, we are considering all decay products to be 
sterile). 

In the absence of decay or oscillations, the recoil elec- 
tron spectrum from neutrino-electron scattering in SK 
can be computed by convolving the 8 B neutrino spec- 
trum pf| with the v e — e elastic scattering cross section 
da/dT (see, e.g., Ref |J), and smearing with the SK 
energy resolution |22]] . For a given decay or oscillation 
scenario, the spectrum is calculated in the same way, 
but now allowing an energy-dependent survival proba- 
bility for the neutrinos. The contribution (about 6 times 
smaller, and with a slightly different differential cross sec- 
tion) from v^ — e elastic scattering must also be included. 
These two spectra, binned every 0.5 MeV in the recoil 
electron total energy E e = T e + m e can then be divided 
bin-by-bin. The resulting "ratio spectrum" Rsk closely 
approximates how SK presents their data, though they of 
course use a much more detailed Monte Carlo to model 
the detector response. Also, in their case the numerator 
in the ratio is the measured counts per bin. 



A. Simple V2 Disappearance 

In the upper panel of Fig. ^| we display the v e survival 
probability as a function of neutrino energy, for a range 
of decay rates. We have chosen parameters in the LMA 
allowed region. The solid curve corresponds to the case of 
no decay. In the lower panel, we show the corresponding 
SK electron energy spectrum, which also includes the 
component. At large energies, the height of the solid lines 
is set just by sin 2 26 = 0.9; SNO expects 0.34 and SK 
expects 0.34 + 0.66/6 = 0.45. The upturn at low energies 
is set by 5m 2 — 4 x 10~ 5 eV 2 . Note that the chlorine 
experiment is mostly sensitive to the 8 B neutrinos, so 
also expects ~ 1/3. In order to correctly calculate the 
rate in the gallium detectors, one must also integrate 
over the finite source region in the Sun. One can see 
from the figure that as the neutrino lifetime approaches 
10~ 4 s/eV or shorter the spectral shape begins to display 
a significant deviation from the flat SK ratio spectrum. 
The total flux also begins to depart significantly from 
the measured value, though only at about the 2a level, 
taking into account the 8 B flux uncertainty. 
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FIG. 2: In the upper panel, the electron neutrino sur- 
vival probability is shown versus neutrino energy. In the 
lower panel, the ratio of measured to expected spectra in 
SK is shown versus the recoil electron total energy, as de- 
tected in neutrino-electron scattering above 5 MeV. The solid 
lines correspond to the LMA solution (5m 2 — 4 x 10 -5 eV 2 , 
sin 2 28 = 0.9) with stable V2- The dashed lines, in order 
of decreasing height, correspond t o z/g lifetimes of r/m = 
10" 3 , 10" 4 ,and 10" 5 s/eV (by Eq. |T|), these correspond to 
decays at typical energies of 0.5, 5, and 50 MeV, respectively). 
The 1258-day data from SK are shown in the lower panel, as 
is the size of the flux normalization uncertainty. The decay 
products of V2 are considered to be sterile. 



B. Compensating Effects from Oscillations 

In Fig. H we used the flatness of the SK spectrum to es- 
timate a bound on the decay rate of v^. However, things 
may not be as straightforward as this basic case. In par- 
ticular, since oscillation effects may lead to an energy 
distortion, we must address the possibility that a cancel- 
lation between oscillation and decay survival probabili- 
ties may produce an acceptably flat spectrum. 

For example, we could choose a 5m 2 somewhat larger 
than those in the allowed LMA region, so that the v e sur- 
vival probability begins to rise toward the lower energy 
end of the SK energy range. Since the decay rate is larger 
for lower energy neutrinos, the question then becomes 
whether the decay might be "just so" in the sense that 
it conspires to cancel this tilt in the survival probability, 
resulting in a spectrum which is apparently flat. Such a 
scenario has recently been addressed in Ref . [||[ |j| . 
We view this scenario as somewhat unnatural, however 
we agree that it cannot be excluded. 

In Fig. [|, we show an example of such a conspiracy, 
using oscillation parameters for which, in the absence of 
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FIG. 3: An example of the parameter degeneracy between 
the lifetime and the mixing parameters is illustrated. The 
description is as for Fig. ^, except that now Sm 2 = 12 x 10~ 5 
eV 2 , which flattens Rsk for r/m = 10 -4 s/eV, making it con- 
sistent with the SK spectral shape. It is about la discrepant 
in terms of normalization. 



decay, the SK spectrum is not very flat. Here, if we 
choose a lifetime of order 10 -4 s/eV we see that a rea- 
sonably flat spectrum may be regained. Note however, 
the large difference between the solid and dashed curves 
in this figure, both in terms of the spectral shape and 
the overall rates. It is important to realize that Kam- 
LAND will essentially break the degeneracy between the 
oscillation parameters and the neutrino lifetime, giving a 
prediction for the mixing effects (without decay) on the 
solar neutrino spectrum. One can then go and look for a 
deviation from this prediction in the solar neutrino spec- 
trum. As is clear from Figs. || and |j| a lifetime of order 
10~ 4 s/eV defines the scale at which decay could possibly 
be distinguished. 



C. Appearance of Active Daughter v\ 

Another possible complication is the possible detec- 
tion of neutrinos produced as decay products. 2 The re- 
placement of a parent neutrino with an active daughter 
may partially hide the effects of decay. Unlike the sit- 



2 The possible detection of decay products was actually noted 
in the first paper on neutrino decay, Ref. H|. However, active 
daughter neutrinos (as opposed to antineutrinos) have since re- 
ceived very little attention in the literature. 



uation discussed above, this would not be sensitive to 
a conspiracy between parameters. Rather, the appear- 
ance and possible detection of active decay products is a 
quite generic expectation, as most plausible decay mod- 
els will feature a neutrino or antineutrino in the final 
state. There are a range of model-dependent possibilities 
for these final-state neutrinos, such as whether the neu- 
trinos are Dirac or Majorana particles and whether the 
decay products are active or sterile [ p6| . In this subsec- 
tion we study the effects of active daughter neutrinos. 



If active daughters are produced, the overall neutrino 
mass scale is important. In virtually all studies of neu- 
trino decay, it has been assumed that the neutrino mass 
spectrum is hierarchical, and thus the mass of the final 
state neutrino may be neglected. This is important, since 
the energy difference between initial and final state neu- 
trinos depends on the values of the masses. If a neutrino 
decays to a massless particle and another neutrino, the 
fraction of the initial neutrino energy carried by the final 
state neutrino ranges from a maximum of Ef /Ei ~ 1 to 
a minimum of Ef/Ei ~ m\/m\. Where a hierarchical 
mass spectrum is a good approximation, the two decay 
products will each carry a large fraction of the original 
neutrino energy, hence the daughter neutrino will, on av- 
erage, be degraded in energy. However, if the neutrino 
masses are degenerate, the daughter neutrino must have 
approximately the same energy as the parent neutrino. 
For example, if we take Sm 2 ~ 4 x 1CP 5 eV 2 and the 
degenerate mass larger than about mi ~ m,2 ~ 0.03 eV 
the energy degradation of the daughter neutrino would 
be less than about 0.5 MeV - the width of the SK energy 
bins. 



Given this, one may wonder if it is possible to set any 
decay bound at all. For example, if the mixing angle 
were exactly maximal, any decay of i>2 to v\ in which the 
neutrino energy was not degraded would be undetectable. 
However, the solar v e survival probability appears to be 
less than 1/2 in the SNO/SK and CI energy region (and 
somewhat larger in the Ga region) strongly disfavoring 
exact maximal mixing. As long as the mixing angle is 
not exactly maximal, decay will cause the relative fluxes 
of v e and to vary as a function of energy, which would 
be distinguishable in both SK and SNO. 

In Fig. |i| we demonstrate the effect on the spectrum if 
Vi decays to the orthogonal state y±, in the limit where 
the daughter neutrinos carry the full energy of the par- 



ent. The expressions in Eq. (3.4) have now to be re- 



3 The possibility that the entire flux (across the whole energy 
range) has decayed to the ui state by the time the neutrinos 
reach Earth would give an energy-independent suppression of 
the solar u e flux, contrary to that observed. 
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FIG. 4: The effect of vi decay to active v\ daughters is 
shown. We have assumed mi and m,2 are nearly degenerate, 
so that the daughter energy is approximately the full energy 
of the parent. The description is as for Fig. 0, except that now 
the dotted lines, in order of increasing height, correspond to 
V2 lifetimes of r/m = 10" 3 , 10~ 4 , and 10~ 5 s/eV. Again, the 
case r/m = 10~ 4 s/eV is discrepant in terms of shape and 
normalization. 
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in order to include the v\ produced by the decay of v%. 
Note that there is no quantum-mechanical interference 
between v\ originally in the beam with those produced 
by the decays of i>2 WA - I n addition, we may take the 
daughter neutrinos to be collinear with the parent neu- 
trinos JUj. 

An interesting feature of Fig. [| is the direction of the 
deviation from a flat spectrum. Rather than a deviation 
downward (as observed in Fig. |J) caused by the depletion 
of the flux at lower energies, we instead obtain a deviation 
upward. This is due to the replacement of a portion of 
the V2 flux by v\ , which has a larger v e component (since 
Qsoiar — 36°). Note, however, that the size of r/m where 
decay effects show up as a significant deviation (either to 
the flat spectrum, or to the total rates) is comparable to 
that shown in Fig. ||. Therefore, the lifetime limits for the 
cases of decay to steriles and decay to active neutrinos of 
non-degraded energy could be comparable (if m\ ~ m,2 
were known). 
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FIG. 5: The effect of U2 decay to active v\ daughters is shown 
again, this time with a "dark-side" angle of 54° instead of 
36° (but still with 8m 2 = 4 x 10" 5 eV 2 , sin 2 26 = 0.9). The 
description is as for Fig. ^, except that now the dotted lines, 
in order of decreasing height, correspond to V2 lifetimes of 
r/m = 10 -3 , 10 -4 , and 10 -5 s/eV. In this unusual case, a 
lifetime as short as 10~ 5 s/eV is required to undo the 1/2/fi 
reversal in v e caused by choosing a "dark-side" angle. 



D. Dark-Side Inversion Undone by Decay 



We now consider an unusual case that requires decay. 
Suppose the solar neutrino parameters live on the so- 
called "dark-side" |28| of the parameter space, that is, 
the hierarchy is inverted such a way that Vi has a larger 
v e component than does v\ . This situation does not pro- 
vide a good fit to the solar neutrino data, as an MSW 
resonance would not take place in the Sun and it is not 
possible to obtain a v e survival probability less than one 
half. However, if we add fast decay to this dark-side 
solution, we can convert the entire flux to v\, obtain- 
ing a solution that is in good agreement with the SK, 
SNO and CI data. That is, decay effectively undoes the 
i>2 — v\ reversal in v e . We have plotted such an example 
in Fig. [|. Note that the curve which most closely resem- 
bles the measured SK data is the one which corresponds 
to the largest decay rate. The only problem with this 
solution lies with Ga flux measurement of greater than 
a half - a difficulty which, alone, is probably sufficient 
to rule out such a solution. So, although such a scenario 
seems unlikely, it is enough, however, to give one pause 
for thought. 



E. Other Possibilities 

In Fig. H we have considered decay products that are 
completely undetectable (this could also be achieved if 
the daughter neutrino energies were artificially taken to 
zero), while in Fig. |] we considered active daughters of 
the full energy, such that their detectability is optimized. 
It is clear that the case where the daughter neutrino 
should be active but have degraded energy, lies some- 
where between the examples considered in Fig. |^ and 
Fig. ^. This case would occur if ra\ <C Tri2- However, 
note that this directly modifies the absolute spectrum, 
which is steeply falling (see Ref 0); the ratio spectra 
shown above cannot be averaged by eye. By comparing 
Figs. H and ||, it should be clear that it could be quite 
difficult to discern such a decay using either the spectral 
shape or the total flux 

Another example is that of the inverted hierarchy, 
where the solar neutrino flux can decay to the eigenstate 
that consists mainly of and v T . Since these flavors 
have a cross-section approximately 6 times smaller than 
v e , this situation is closer to the sterile daughters exam- 
ple presented in Fig. H than to the active daughter v\ 
case presented in Fig. HI 

The decay products might include antineutrinos rather 
than neutrinos or, more generally, a mixture of the two. 
Since v ^ and V T have similar neutrino-electron scattering 
cross-sections to and v T , they would also be difficult 
to detect. 

In Ref. 0, p9| , the appearance of v e as a decay prod- 
uct is discussed. An estimated limit on the solar v e flux 
has been obtained for the SK data in Ref. |30| , using the 
technique of Ref. fl3l]] , and is quoted as 3.5% of the 8 B 
v e flux (at the 95% C.L.). Ref. Q takes this to imply a 
stringent decay limit of r/m ~ 10 -3 s/eV. Though they 
state that the v\ carries less energy than the parent neu- 
trino, it appears that the bound is derived without taking 
this into account. Indeed, as we have stressed above, the 
daughter may carry nearly the full energy of the parent. 
But if not, it makes a big difference for the yield, since 
the v e p — > e + n cross section is nearly quadratic in neu- 
trino energy. And the neutrino spectrum is falling steeply 
with energy, so decay products from a high energy can 
be hidden at a lower energy. Also, if the daughter energy 
is less than the parent energy, then the limit as quoted 
from Ref. |3(| does not apply. The limit of 3.5% assumes 
that the v e spectrum has the same shape as the 8 B v e 
spectrum. If that assumption is relaxed, then the flux 
limit is more conservatively about 10%. Thus the decay 
limit is certainly considerably weaker than r/m ~ 10~ 3 
s/eV. 

Finally, one might reasonably ask if the solar neutrino 
data can be explained by neutrino decay. The formula- 
tion of this question that seems to be the most interesting 
is this: Can the solar v e and LSND |32| v e data be made 
consistent? The combined solar, atmospheric, and LSND 
results require three independent Sm 2 values, whereas 
only two are allowed in three-generation mixing. Thus, 
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our question can be rephrased: Can one explain all of the 
data with the two Sm 2 values and one r/m value? Let us 
assume the LSND vacuum mixing parameters. For the 
large 8m 2 ~ 1 eV 2 , matter effects in the Sun are negligi- 
ble. We are free to invert the sign of dm 2 so that the solar 
v e ~ V2 (the LSND mixing angle is very small). Suppose 
decay turns this into ~ 1/3 v% + 2/3 v\ ~ 1/3 v e + 2/3 
in the SK energy region, roughly matching the SNO and 
SK observations (we ignore the spectral distortion) . But 
then decay is complete in the gallium-detector energy 
range, predicting nearly zero v e flux there, in gross dis- 
agreement with observations. 



IV. OTHER NEUTRINO DECAYS 

Let us suppose that a model-independent limit on so- 
lar neutrino decay can be arrived at, and that it reaches 
the scale of r/m < 10~ 4 s/eV. Once such a limit has 
been established, can it be used to set meaningful limits 
on the possible decay of neutrinos from other sources? 
In particular, we shall consider the decay of atmospheric 
neutrinos. It is immediately obvious that if we restrict 
ourselves to the three active neutrinos, it is difficult to 
arrange for decay to play a role in the atmospheric neu- 
trino problem. Since two of the three mass eigenstates 
have large v e components, it is hard to see how decay 
could alter the atmospheric flux, while leaving the v e 
flux unaltered, as the data suggests. 

We can make this point more precise. Let us define 
the neutrino mass eigenstates to be such that m 3 > 
TO2 > mi. Considered as a function of mass squared, 
in a normal hierarchy two states are close together (the 
solar 5m 2 ), and are well below the third state (by the at- 
mospheric 8m 2 ). The sign of the solar Sm 2 is fixed by the 
observation of matter effects in the Sun (for the opposite 
sign, there is no MSW resonance). However, the sign of 
the atmospheric 8m 2 that dictates vacuum oscillations is 
unknown. An inverted hierarchy is obtained when this 
state is well below the solar doublet on a mass-squared 
scale. See also Fig. but note that the scale there is 
mass, not mass squared. 

If we have a normal hierarchy, the solar decay limit 
applies to T21/TO2, whereas if we have an inverted hierar- 
chy, the limit applies to both r 3 2/m 3 and r 3 i/m 3 . In the 
inverted case, only v 2 — > v\ is not directly constrained 
(in the inverted case, we relabel so that m 3 > m 2 > m\). 
However, since this decay mode has virtually the same 
phase space as z^ 3 — > u\, it must also be constrained, 
unless there were a large hierarchy in the respective cou- 
plings. It is important to remember that the flavor con- 
tent of the mass eigenstates is irrelevant for their decays. 

In the case of a normal hierarchy, we can translate the 
bound on v 2 decay into one on z/ 3 if we make certain 
assumptions. Since the quantity limited by experiment 
is t jm 2 and r/m 2 ~ (g 2 Sm 2 )~ 1 (using Eq. (pH|)), then 
the overall mass scale is mostly irrelevant (at least for 
the hierarchical case) for setting a bound on j. A solar 



decay limit of 10 4 s/eV would translate into a limit on 
the coupling of 



S21 < 3 x 10- 



/10- 5 eV 2 
\ 8m 2 



(4.1) 



For the solar LMA solution, Sm 2 ~ 4 x 10~ 5 eV 2 , one 
would obtain a limit of 921 < 10~ 5 , slightly more re- 
strictive than the limit from meson decays of g 2 < 10~ 4 . 
However, from Section II we see that the correspondence 
between neutrino lifetime and the coupling constant de- 
pends sensitively on the mass hierarchy. 

In order to translate a bound on v 2 into one on ^ 3 , 
we would have to assume that g is approximately univer- 
sal among generations. Of course, this is a very model- 
dependent assumption. The limit on the v 3 lifetime 
would become 



TJL ^ 321 x Srr& 1Q _ 6 
m 3 m 2 tajj 



(4.2) 



The natural scale for setting a limit on neutrino de- 
cay with atmospheric neutrinos is r/m ~ L/E ~ 
10 4 km/300 MeV ~ 10~ 10 s/eV. Thus if there were sig- 
nificant decay in the atmospheric neutrinos, a significant 
hierarchy in the couplings ((731/321) 2 ~ 10 4 would be re- 
quired. The exception to this is to greatly increase the 
Sm 2 , as in Ref. |3j|. This model is now strongly dis- 
favored since it cannot accommodate large-angle mixing 
with an active neutrino in the solar sector. Finally, as 
noted, if there is an inverted hierarchy, a limit on solar 
neutrino decay may apply directly to atmospheric neu- 
trino decay, since the same states may be involved. 

A limit on solar neutrino decay would also apply to 
decays to "phantom" neutrinos. Light sterile neutrinos 
can be used to add new mass eigenstates anywhere rel- 
ative to the standard hierarchy. These new mass eigen- 
states might be inaccessible by flavor mixing (or with very 
small angles), but might be reached via decays (which 
connect mass eigenstates directly). If so, the presence of 
these phantom neutrinos could result in an apparent non- 
unitarity of the 3x3 mixing matrix. Indeed, such tests 
could infer their existence even if decay were not seen 
directly. These phantom neutrinos can also be lower in 
mass than v\ in the standard case; however, the very long 
SN1987a lifetime limit |34 would apply to those states. 



V. CONCLUSIONS 

Solar and supernova neutrinos have been observed on 
Earth suggesting, naively, that they do not decay. How- 
ever, this conclusion cannot immediately be drawn unless 
one can rule out certain subtleties. When considering 
these possibilities we must be careful to distinguish the 
mass eigenstates (for which the lifetimes are properly de- 
fined) and the flavor eigenstates. 

1. While the LMA solution is an excellent fit to the 
solar neutrino data, one cannot completely rule out 
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more exotic scenarios until KamLAND con- 
firms the LMA parameters in an experiment that 
does not rely upon properties of the Sun. Kam- 
LAND will determine the "solar" mixing parame- 
ters using antineutrinos rather than neutrinos, in 
vacuum rather than in matter, and using a much 
shorter baseline than the Earth-Sun distance. In 
addition, exotic effects such as flavor-changing neu- 
tral currents, or resonant spin-flip transitions can 
be eliminated. Only then shall we have the com- 
plete confidence in the LMA solution necessary to 
fully utilize the potential of the solar neutrino beam 
as a probe of non-standard neutrino properties. 

2. The flat energy spectrum observed in SK and SNO 
is well described by the LMA solution, which would 
seem to argue against an energy-dependent distor- 
tion of the survival probability, as would be charac- 
teristic of decay. However, for parameters outside 
the LMA region one can obtain an energy depen- 
dent oscillation survival probability to offset the ef- 
fects of decay. This possibility may be eliminated 
by KamLAND. 

3. Decay will typically (though not always) cause a 
depletion of the total solar neutrino flux. While this 
is also a possible way to identify decay, we should 
keep in mind the 20% uncertainly in the 8 B flux 
normalization. That uncertainty will be reduced 
by future neutral-current data from SNO, as well 
as direct nuclear-physics measurements |36|. 

4. In the limit that the neutrino masses are degen- 
erate, a daughter neutrino produced by decay will 
carry the full energy of the parent neutrino, and 
could be detected in solar neutrino experiments. 
The replacement of the parent neutrino with an 
active daughter of the same energy could obscure 
the characteristic features of decay. This is espe- 
cially pertinent in the case of the decay V2 — > vi , 



where both v-i and v\ have large v e projections. 

5. Decays producing V e daughters should be readily 
detectable, provided the energy of the V e is not too 
degraded. If the hierarchy is inverted, there may 
be decays to the state that is dominantly and 
v T (or the antiparticle state). These are harder to 
detect in SK, though they would add to the integral 
neutral-current rate in SNO. 

6. If a model-independent limit on neutrino decay can 
be established, the bound will likely be of order 
r/m > 10~ 4 s/eV. This is close to the limit ob- 
tained via meson-decay bounds (which may or may 
not apply; see above). Though the bound is quite 
weak, the solar neutrino flux would, however, still 
provide the best limit on non-radiative neutrino de- 
cay. This limit is about eight orders of magnitude 
too weak to rule out the decay of astrophysical neu- 
trinos (e.g. , from a Galactic supernova or a distant 
AGN pTfl) on their journey to Earth. We thus have 
absolutely no guarantee we shall be able to use such 
neutrinos to probe the astrophysics of these sources 
without taking decay into account. 
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